The nanofluidics can explain 
ascent of water in tallest trees 

Henri GouirQ 
Abstract. 

In Amazing numbers in biology, Flindt reports a giant, 
128 meter-tall eucalyptus, and a 135 meter-tall se- 
quoia 1 . However, the explanation of the maximum 
altitude of the crude sap ascent and consequently 
the main reason of the maximum size that trees can 
reach is not well understood 2 . 

According to tree species, the crude sap is driven in 
xylem microtubes with diameters ranging between 
50 and 400 micrometers. The sap contains diluted 
salts but its physical properties are roughly those of 
water; consequently, hydrodynamic, capillarity and 
osmotic pressure yield a crude sap ascent of a few 
tens of meters only 3 . 

Today, we can propound a new understanding of the 
ascent of sap to the top of very tall trees thanks to 
a new comparison between experiments associated 
with the cohesion-tension theory 4 and the disjoining 
pressure concept 5 . 

Here we show that the pressure in the water-storing 
tracheids of leaves can be strongly negative whereas 
the pressure in the xylem microtubes of stems may 
remain positive when, at high level, inhomogeneous 
liquid nanolayers wet the xylem walls of microtubes. 
The nanofluidic model of crude sap in tall trees 6 
discloses a stable sap layer up to an altitude where 
the pancake layer thickness 7 coexists with the dry 
xylem wall and corresponds to the maximum size of 
tallest trees. In very thin layers, sap flows are widely 
more significant than those obtained with classical 
Navier-Stokes models and consequently are able to 
refill stomatic cells when phloem embolisms super- 
vene. 

These results drop an inkling that the disjoining 
pressure is an efficient tool to study biological liquids 
in contact with substrates at a nanoscale range. 



The most classical explanation of the sap ascent phe- 
nomenon in tall trees is given by the cohesion-tension 
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theory propounded in 1894 by Dixon and Jolj3, fol- 
lowed by a quantitative analysis of the sap motion pro- 
posed by van der Honert in 194P: According to this 
theory, the crude sap fills water-tight microtubes of 
dead xylem cells and its transport is due to a gradient 
of negative pressure producing the traction necessary to 
lift water against gravity. A main experimental check- 
ing on the cohesion-tension theory comes from an appa- 
ratus called Scholander pressure chamber (see Fig.[T]). 
The decrease in the negative pressure is related to the 
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Figure 1: Sketch of the Scholander pressure cham- 
ber (or Scholander pressure bombj^l . A leaf attached to a 
stem is placed inside a sealed chamber. Compressed air is 
slowly added to the chamber. As the pressure increases to 
a convenient level, the sap is forced out of the xylem and 
is visible at the cut end of the stem. The required pressure 
is opposite and of equal magnitude to the water pressure in 
the water-storing tracheids in the leaf. 

closing of the aperture of microscopic stomata in leaves 
through which water vapour is lost by transpiration. 

Nonetheless, several objections question the 
cohesion-tension theory: 

We first refer to the well-known book by M.H. Zim- 
mermanrP. He said: "The heartwood is referred to 
as a wet wood, ft may contain liquid under positive 
pressure while in the sapwood the transpiration stream 
moves along a gradient of negative pressures. Why is 
the water of the central wet core not drawn into the 
sapwood? Free water, i.e. water in tracheids, decreases 
in successively older layers of wood as the number 
of embolized tracheids increases. The heartwood is 
relatively dry i.e. most tracheids are embolized. It is 
rather ironic that a wound in the wet wood area, which 
bleeds liquid for a long period of time, thus appears to 
have the transpiration stream as a source of water, in 
spite of the fact that the pressure of the transpiration 
stream is negative most of the time! It should be quite 
clear by now that a drop in xylem pressure below 
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a critical level causes cavitations and normally puts 
the xylem out of function permanently. The cause of 
such a pressure drop can be either a failing of water 
to the xylem by the roots, or excessive demand by 
transpiration." 

At great elevation in trees, the value of the negative 
pressure increases the risk of cavitation and conse- 
quently, the formation of embolisms may cause a 
breakdown of the continuous column of sap inducing 
the leaf death. The crude sap is a liquid bulk with a 
superficial tension a lower than the superficial tension 
of pure water which is about 72.5 cgs at 20° Celsius; 
if we consider a microscopic gas- vapour bubble with a 
diameter 2R smaller than xylem microtube diameters, 
the difference between the gas-vapour pressure and the 
liquid sap pressure is expressed by the Laplace formula 
P v — Pi — 2 a/7?; the vapour-gas pressure is positive 
and consequently unstable bubbles must appear when 
R > —2a/ Pi. For a negative pressure Pi = —0.6 MPa 
in the sap, we obtain R > 0.24 /tm; then, when all the 
vessels are tight filled, germs naturally pre-existing in 
crude water may spontaneously embolize the tracheids. 
HaberlandP 1 described water-storing tracheids in 
leaves; they are roundish in shape, and located either 
at the tips of the veins or detached from transporting 
xylem. In more recent papers they have been called 
tracheid idioblat^^. The spacing considered in Prid- 
geon's paper—' is about 2 /im or less at the top of trees 
as suggested on Fig. 2 of the paper by Koch et oP and 
has a good size to prevent cavitation for nucleus germs 
of the same order of magnitude. 

No vessels are continuous from root to stem, from stem 
to shoot, and from shoot to petiole. The vessels do not 
all run neatly parallel and form a network generally up 
a few centimeters long. The ends usually taper out; 
it is very important for the understanding of water 
conduction to realize that the water does not leave 
a vessel in axial direction through the very end but 
laterally along a relatively long stretch where the two 
vessels, the ending and the continuing ones, run side 
by side. 

The vascular bundles of some leaves are surrounded 
by a bundle sheath, containing a suberized layer 
comparable to the one of the Casperian strip in the 
roots^. This seal separates the apoplast into two 
compartments, one inside and the other one outside 
the bundle sheath. The two areas are only connected 
by the plasmodesmata that connect living cells. The 
pressure in the intact, water-containing neighbouring 
tracheids, may still be negative; a considerable pres- 
sure drop therefore exists across the pit membranes. 
Pressure chamber measurements cannot be considered 
as pressure values of the stem xylem without special 
precautions, simply because they are taken elsewhere. 



Hydraulically then, the leaf is very sharply separated 
from the stem. The wet wood area of elms appears to 
act like a single, giant osmotic cell that is separated 
from the sapwood area by a semi permeable mem- 
brane. This can be visualized somewhat like a Traube 
membrane, as early plant physiologists called fp^. 
An other objection in the perfect confidence to the 
cohesion-tension theory was the experiments by Pre- 
ston^ who demonstrated that tall trees survived by 
overlapping double saw-cuts made through the cross- 
sectional area of the trunk to sever all xylem elements. 
This result, confirmed by several authors (e.gp^), does 
not seem in agreement with the possibility of strong 
negative pressures in the water-tight microtubes. 
Using a xylem pressure probe, Balling et ar-^ showed 
that, in many circumstances, this apparatus does not 
measure any water tensiorP^. 

As a consequence of these various pieces of evidence, 
the main question is: 

The negative pressure measured by the Sholander pres- 
sure chamber being the pressure in the water-storing 
tracheids, is it possible that the pressure in the xylem 
microtubes of the stem remains positive? 

A positive answer to this question comes from the 
concept of disjoining pressure which is able to inter- 
pret very thin vertical films of liquid wetting solid sub- 
strates. The simplest review relative to the disjoining 
pressure is presented in the well-known monograph by 
Derjaguin et aW. An apparatus measuring the disjoin- 
ing pressure denoted II is due to Sheludko^. The ap- 
paratus is schematically described on Fig. [2] and allows 
to understand the disjoining pressure concept. Fluids 
and solids are at the same temperature. The film is 
thin enough for the gravity effect across the layer to 
be neglected. The hydrostatic pressure in a thin liquid 
layer included between a solid substrate and a vapour- 
gas bulk differs from the pressure in the liquid bulk 
contained in the reservoir. The forces arising from the 
thinning of a film of uniform thickness h yield the dis- 
joining pressure as a function of h (II = 11(h)). The 
disjoining pressure is equal to the difference between 
the pressure P Vb through the interfacial surface of the 
thin liquid layer (which is the total pressure of gases 
and vapour bulk) and the pressure Pi b at the top of 
the liquid bulk with dissolved gases contained in the 
reservoir from which the thin liquid layer extends: 

n = Pv b -Pt„. 

As a vapour-gas pressure, the pressure P Vb is always 
positive; depending on H-values, the pressure Pi b may 
be negative. 

These comments illuminate an original comparison 
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Figure 2: The Sheludko apparatus . Schematic dia- 
gram of the technique to determine the disjoining pressure 
isotherms of wetting films on a solid substrate: a circular 
wetting film is formed on a flat substrate to which a microp- 
orous filter is clamped. A pipe connects the filter filled with 
the liquid to a reservoir containing the liquid bulk that can 
be moved by a micrometric device. The disjoining pressure 
is equal to il = (p lb — p v ) g H, where g is the acceleration 
of gravity and p lb , p Vb are the densities in the liquid and the 
vapour-gas bulks, respectively (Fronfflpage 332). 



between the Sholander pressure chamber experiment 
presented in Fig. [T] and the apparatus proposed by She- 
ludko in Fig. [5] and make a total change in the inter- 
pretation of Sholander pressure chamber data for the 
tallest trees: 

In the Sholander pressure chamber, the cohesion- 
tension theory assumes that, at a given level, liquids 
entirely convey the pressure as is the case for incom- 
pressible fluids. This assumption is for thin layers at 
odds with the Sheludko experiment where the disjoining 
pressure highlights a strong difference between liquid 
bulk and thin layer pressures. The xylem of stem walls 
is associated with the solid substrate and the water- 
storing tracheids in leaves are associated with the liquid 
bulk reservoir of Sheludko's experiment. 
In our new interpretation, the negative pressure mea- 
sured in the Sholander pressure chamber is the pressure 
in water-storing tracheids corresponding to the liquid 
bulk reservoir of Sheludko while the pressure in the 
xylem microtubes of stems remains positive as in the 
wetting thin layer of Sheludko's experiment when H- 
value is equal to the level in the tree. 

As proven irP, Chapter 2, the Gibbs free energy per 
unit area of the liquid layer denoted G can be expressed 
as a function of h thanks to the relation: 

r+oo 

G{h) = / U(u)du, 

Jh 

where h = is associated with a dry wall in contact 
with the vapour-gas bulk and h = +oo is associated 
with a wall in contact with the liquid bulk when the 
value of G is equal to 0. The coexistence of two film seg- 



Figure 3: The pancake layer thickness . The con- 
struction of the tangent to curve G(h) issued from point W 
of coordinates (0, G(0)) yields point P; point W is associ- 
ated with a high-energy surface of the dry wall and point P 
is associated with the pancake thickness hjpl . 



ments with different thicknesses is a phenomenon which 
can be represented by the equality of chemical potential 
and superficial tension of the two films^. A spectacular 
case corresponds to the coexistence of a liquid film of 
thickness h p and the dry solid wall; the liquid film is the 
so-called pancake layer and corresponds to the minimal 
thickness for which a stable wetting film damps a solid 
walPa 

This minimal thickness h p verifies: 



G(0) = G(h p ) + MI(/i p ). 



(1) 



The geometrical interpretation of Eq. |T]) is proposed 
on Fig.|3 

Now, we consider a film of thickness h x at level x. Only 
liquid water films of thicknesses h x > h p are stable. The 
disjoining pressure of the mixture of water and perfect 
gas i s the same as for a single liquid far under its critical 
pohrP^. Calculations limited to overstrained liquid 
layers of a few nanometers thick yield the disjoining 
pressure in the approximated for m 20 * 21 !: 



U(h x ) 

2 A 
pi 

[(7 2 + 7 3 )7 2 Pi + (7i - 7 2 ft)(7 2 - 7 3 ) e " hxT ] 



-~ [(7i ~ 7 2 ft)(7 2 + 7 3 ) e ' , ' a!T + (72 ~ 7 3 )7 2 ftl 



[(7 2 +l 3 ) 2 e h * r - (7 2 -7 3 ) 2 e- h * T ] 2 



(2) 



where p { is the liquid density in normal conditions, q is 
the isothermal sound velocity in the liquid water bulk, 
7i> 7 2 i 7 3 an d T are posit ive co nstants given by the 
mean field molecular theorjP^^: 
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Constants cu and Ci a are associated with Hamaker co- 
efficients of interaction of liquid VS liquid and liquid vs 
solicP; a; denotes the fluid molecular diameter and 
S = i(ai + a 3 ), where a B denotes the solid molecular 
diameter; mi, m s denote the masses of fluid and solid 
molecules; p aol denotes the solid density. Expression 
@ differs from Lifshitz's relation^ where the disjoining 
pressure of microscopic layers of liquid, assumed to be 
incompressible, has a behavior in the form H(h) ~ h~ 3 . 
To obtain the pancake thickness corresponding to the 
smallest thickness of the liquid layer, we draw the 
graphs of G(h x ) and H(h x ) when h x £ [(1/2) crt,£], 
where t is a length of a few tens of angstroms. Let us 
note that d = 1/r is the natural reference length scale. 
For a few nanometers, the film thickness is not exactly 
h x ; we must add the thickness estimated at 2 <t; of the 
liquid part of the liquid-vapour interface bordering the 
liquid layer and the nanolayer thickness is approxima- 
tively e x « h x +2aF%. 

Our aim is now to point out a numerical example 
such that previous results provide a value of maxi- 
mum height for a vertical water film wetting a wall 
of xylem. We considered water at 20° Celsius. The 
experimental estimates of coefficients are expressed in 
c.g.s. unitsP 1 ^ 1 : p l = 0.998, c ; = 1.478 x 10 5 , 
cu = 1.4 x 10" 58 , oi = 2.8 x 10" 8 , mi = 2.99 x 10" 23 . 
We deduce 7 2 = 54.2, 7 3 = 506, d = 2.31 x 10" 8 . 
We consider the Young contact angle between the 
xylem wall and the liquid-vapour water interface as 
« 50° (this value is an arithmetic average of different 
Young angles proposed in the literature^). Coefficients 
cis and Yj can be obtained from the substrate-liquid 
surface free energy expressed in the forrrP"^: 



0(P S 



1 2 
-llPs + o 7 2 P 3 - 



Here p 3 ~ p l denotes the fluid density value at the sur- 
face; from the superficial tension a and Young's condi- 
tion, we immediately get y 1 as 75 (p^). 
In the upper graph of Fig. [4] we present the free energy 
graph G(h x ). Due to h x > (1/2) ai, it is not numeri- 
cally possible to obtain the limit point W corresponding 
to the dry wall; point W is obtained by an interpolation 
associated with the concave part of the G-curve. Point 
P follows from the drawing of the tangent line issued 
from W to the G-curve. In the lower graph of Fig.|4]we 
present the disjoining pressure graph Tl(h x ). The phys- 
ical part of the disjoining pressure graph corresponding 
to dH/dh x < is associated with a liquid layer of sev- 
eral molecules thick. The part d\l/dh x > is also 
obtained by Derjaguin et The reference length d 
is of the same order as ai and is a good length unit for 
very thin films. The total pancake thickness e p is of 
one nanometer order corresponding to a good thickness 
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Figure 4: The maximum height of trees . Upper 
graph: G(h x )-gra,ph. The unit of x— axis graduated by h x 
is d = 2.31 X 10 — 8 cm ; the unit of y— axis is one c.g.s. unit 
of surface tension. Lower graph: n(h^)-graph. The unit of 
x— axis graduated by h x is d = 2.31 X 10 — 8 cm; the unit of 
y— axis is one atmosphere. 



value for a high-energy surface^"^; consequently in the 
tall trees, at high level, the thickness of the layer is of 
a few nanometers. The point P of the lower graph cor- 
responds to the point P of upper graph. Let us note 
that the crude sap is not pure water; its liquid-vapour 
surface tension has a lower value than the surface ten- 
sion of pure water and it is possible to obtain the same 
spreading coefficients with less energetic surfaces. 
When x p corresponds to the altitude of the pancake 
layer, II ~ p^xj^^. To this altitude, we add 20 
meters corresponding to the ascent of sap due to capil- 
larity and osmotic pressure and we obtain on the lower 
graph of Fig. [4] a maximum film height of approxima- 
tively 140 meters (20 + 120 meters) corresponding to 
12 atmospheres, which is of the same level order as the 
topmost trees. 

These results arising from molecular physics require a 
comparison between the behaviours of liquid motions 
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both in microtubes and in nanolayers : 

When the xylem microtubes are tight filled with 
crude sap, the liquid motions are Poiseuille flowP and 
to be efficient the microtube radii should be as wide as 
possible, which is not the case. Due to the liquid in- 
compressibility the flow is very rigid and the pressure 
effects are fully propagated onto the microtube walls. 

When the xylem microtubes are partially filled with 
thin water films, there are qualitative experiments for 
slippage on the walls when the film thickness is of the 
mean free path order^"^; the boundary condition on 
the walls writes: 

u — L s du/dn, 

where u is the liquid velocity, du/dn is the normal 
derivative at the wall and L a is the so-called Navier 
lengtf^^. The Navier length may be as large as a few 
microns^ and we obtained the mean liquid velocity u 
along a thin layer frorrP: 

vu — h x + Ls^j [grad U(h x ) - gi] , (3) 

where v denotes the kinematic viscosity and i the unit 
vertical vector. Consequently the slippage condition 
multiplies the flow rate by a factor of (l + 3L„/h x ). For 
example, if h x — 3nm and L s — 100 nm, which is a 
Navier length of small magnitude with respect to exper- 
iments, the multiplicative factor is 10 2 ; if L s = 7 fim, as 
considered irP3, the multiplicative factor is 10 4 which 
is of the same order as for nanotube observations^. 
Equation Q is mainly realistic at the top of tallest 
trees where the xylem network is strongly ramified; the 
heartwood may contain liquid under positive pressure 
and is connected with the sapwoocP. The flow rate 
can increase or decrease due to the spatial derivative of 
the thickness h x and consequently depends on the local 
disjoining pressure value. The tree's versatility adapts 
the disjoining pressure gradient effects by opening or 
closing the stomatic pits, so that the bulk pressure in 
micropores can be more or less negative and so, the 
transport of water is differently dispatched in the stem 
parts. 

The sap motion is induced by the transpiration across 
micropores located in tree leaves^. It seems natural to 
surmise that the diameters of xylem microtubes might 
result of a competition between evaporation which re- 
duces the flow of sap and the flux of transpiration in mi- 
cropores inducing the motion. It is noticeable that if we 
replace the flat surfaces of the microtubes with wedge 
geometry or corrugated surface, it is much easier to ob- 
tain the complete wetting requirement; thus, plants can 
avoid having very high energy surfaces. Nonetheless, 
they are still internally wet if crude sap flows through 



wedge shaped corrugated pores. The wedge does not 
have to be perfect on the nanometric scale to signifi- 
cantly enhance the amount of liquid flowing at modest 
pressures, the walls being considered as plane surfaces 
endowed with an average surface energy. 

Methods. 



We compare two experiments: 

- The Scholander pressure bomb experiment (1955) 
based on the cohesion-tension theory (1894) * n which 
liquids are considered to be incompressible. 

- The Sheludko experiment (1967) based on the concept 
of disjoining pressure in DLVO theory (1948) that 
highlights a strong difference between liquid bulk and 
thin layer pressures. 

The theoretical results allow us to obtain: 

- The computation of tallest trees' level that fits with 
real facts. 

- The interpretation of the motion in xylem microtubes 
by using the shallow water approximation and the 
slippage on walls at the nanometric scale. 
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